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Preface 

These are general notes on tensor calculus originated from a collection of personal notes 
which I prepared some time ago for my own use and reference when I was studying the 
subject. I decided to put them in the public domain hoping they may be beneficial to some 
students in their effort to learn this subject. Most of these notes were prepared in the 
form of bullet points like tutorials and presentations and hence some of them may be more 
concise than they should be. Moreover, some notes may not be sufficiently thorough or 
general. However this should be understandable considering the level and original purpose 
of these notes and the desire for conciseness. There may also be some minor repetition 
in some places for the purpose of gathering similar items together, or emphasizing key 
points, or having self-contained sections and units. 

These notes, in my view, can be used as a short reference for an introductory course on 
tensor algebra and calculus. 1 assume a basic knowledge of calculus and linear algebra 
with some commonly used mathematical terminology. 1 tried to be as clear as possible and 
to highlight the key issues of the subject at an introductory level in a concise form. 1 hope 
1 have achieved some success in reaching these objectives at least for some of my target 
audience. The present text is supposed to be the first part of a series of documents about 
tensor calculus for gradually increasing levels or tiers. 1 hope 1 will be able to finalize and 
publicize the document for the next level in the near future. 
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1 Notation, Nomenclature and Conventions 

• In the present notes we largely follow certain conventions and general notations; most of 
which are commonly used in the mathematical literature although they may not be univer- 
sally adopted. In the following bullet points we outline these conventions and notations. 
We also give initial definitions of the most basic terms and concepts in tensor calculus; 
more thorough technical definitions will follow, if needed, in the forthcoming sections. 

• Scalars are algebraic objects which are uniquely identified by their magnitude (abso- 
lute value) and sign (±), while vectors are broadly geometric objects which are uniquely 
identified by their magnitude (length) and direction in a presumed underlying space. 

• At this early stage in these notes, we generically define “tensor” as an organized array 
of mathematical objects such as numbers or functions. 

• In generic terms, the rank of a tensor signifies the complexity of its structure. Rank-0 
tensors are called scalars while rank-1 tensors are called vectors. Rank-2 tensors may be 
called dyads although this, in common use, may be restricted to the outer product of two 
vectors and hence is a special case of rank- 2 tensors assuming it meets the requirements 
of a tensor and hence transforms as a tensor. Like rank-2 tensors, rank-3 tensors may 
be called triads. Similar labels, which are much less common in use, may be attached to 
higher rank tensors; however, none will be used in the present notes. More generic names 
for higher rank tensors, such as polyad, are also in use. 

• In these notes we may use “tensor” to mean tensors of all ranks including scalars (rank-0) 
and vectors (rank-1). We may also use it as opposite to scalar and vector (i.e. tensor of 
rank-n where n > 1). In almost all cases, the meaning should be obvious from the context. 

• Non-indexed lower case light face Latin letters (e.g. / and h ) are used for scalars. 

• Non-indexed (lower or upper case) bold face Latin letters (e.g. a and A) are used for 
vectors. The exception to this is the basis vectors where indexed bold face lower or upper 
case symbols are used. However, there should be no confusion or ambiguity about the 
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meaning of any one of these symbols. 

• Non-indexed upper case bold face Latin letters (e.g. A and B) are used for tensors (i.e. 
of rank >1). 

• Indexed light face italic symbols (e.g. a,; and Bj k ) are used to denote tensors of rank > 0 
in their explicit tensor form (index notation). Such symbols may also be used to denote 
the components of these tensors. The meaning is usually transparent and can be identified 
from the context if not explicitly declared. 

• Tensor indices in this document are lower case Latin letters usually taken from the 
middle of the Latin alphabet like (■ i,j , k ). We also use numbered indices like (*i, * 2 , • • • , h) 
when the number of tensor indices is variable. 

• The present notes are largely based on assuming an underlying orthonormal Cartesian 
coordinate system. However, parts of which are based on more general coordinate systems; 
in these cases this is stated explicitly or made clear by the content and context. 

• Mathematical identities and definitions may be denoted by using the symbol ‘=’. How- 
ever, for simplicity we will use in the present notes the equality sign to mark identities 
and mathematical definitions as well as normal equalities. 

• We use 2D, 3D and ?rD for two-, three- and n-dimensional spaces. We also use Eq./Eqs. 
to abbreviate Equation/Equations. 

• Vertical bars are used to symbolize determinants while square brackets are used for 
matrices. 

• All tensors in the present notes are assumed to be real quantities (i.e. have real rather 
than complex components). 

• Partial derivative symbol with a subscript index (e.g. i) is frequently used to denote the 
i th component of the Cartesian gradient operator V: 


di = V,; 


d_ 

dxi 


( 1 ) 
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• A comma preceding a subscript index (e.g. , i) is also used to denote partial differentia- 
tion with respect to the i th spatial coordinate in Cartesian systems, e.g. 


A, 


8A 

dxi 


( 2 ) 


• Partial derivative symbol with a spatial subscript, rather than an index, are used to 
denote partial differentiation with respect to that spatial variable. For instance 


d r = V r 


d_ 

dr 


( 3 ) 


is used for the partial derivative with respect to the radial coordinate in spherical coordi- 
nate systems identified by (r, 6 , 0) spatial variables. 

• Partial derivative symbol with repeated double index is used to denote the Laplacian 
operator: 

da = dA = V 2 = A (4) 

The notation is not affected by using repeated double index other than i (e.g. djj or dkk)- 
The following notations: 

dl d 2 d{& (5) 

are also used in the literature of tensor calculus to symbolize the Laplacian operator. 
However, these notations will not be used in the present notes. 

• We follow the common convention of using a subscript semicolon preceding a subscript 
index (e.g. Aki-i) to symbolize covariant differentiation with respect to the i th coordinate 
(see § 6). The semicolon notation may also be attached to the normal differential operators 
to indicate covariant differentiation (e.g. V ; j or d- % to indicate covariant differentiation with 
respect to the index %). 

• All transformation equations in these notes are assumed continuous and real, and all 
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derivatives are continuous in their domain of variables. 

• Based on the continuity condition of the differentiable quantities, the individual differ- 
ential operators in the mixed partial derivatives are commutative, that is 

d idj = <);(): ( 6 ) 

• A permutation of a set of objects, which are normally numbers like (1,2, ...,n) or 
symbols like (■ i,j,k ), is a particular ordering or arrangement of these objects. An even 
permutation is a permutation resulting from an even number of single-step exchanges 
(also known as transpositions) of neighboring objects starting from a presumed original 
permutation of these objects. Similarly, an odd permutation is a permutation resulting 
from an odd number of such exchanges. It has been shown that when a transformation 
from one permutation to another can be done in different ways, possibly with different 
numbers of exchanges, the parity of all these possible transformations is the same, i.e. all 
even or all odd, and hence there is no ambiguity in characterizing the transformation from 
one permutation to another by the parity alone. 

• We normally use indexed square brackets (e.g. [A] i and [V/].) to denote the i th compo- 
nent of vectors, tensors and operators in their symbolic or vector notation. 

• In general terms, a transformation from an nD space to another nD space is a corre- 
lation that maps a point from the first space (original) to a point in the second space 
(transformed) where each point in the original and transformed spaces is identified by n 
independent variables or coordinates. To distinguish between the two sets of coordinates 
in the two spaces, the coordinates of the points in the transformed space may be notated 
with barred symbols, e.g. ( x 1 ,x 2 , . . . ,x n ) or (aq,x 2 , • • • ,x n ) where the superscripts and 
subscripts are indices, while the coordinates of the points in the original space are notated 
with unbarred symbols, e.g. (x 1 , x 2 , . . . , x n ) or (xi, x 2 , . . . , x n ). Under certain conditions, 
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such a transformation is unique and hence an inverse transformation from the transformed 
to the original space is also defined. Mathematically, each one of the direct and inverse 
transformation can be regarded as a mathematical correlation expressed by a set of equa- 
tions in which each coordinate in one space is considered as a function of the coordinates 
in the other space. Hence the transformations between the two sets of coordinates in 
the two spaces can by expressed mathematically by the following two sets of independent 
relations: 

x l — ^(x 1 , x 2 , , x n ) & x 1 = x\x 1 ,x 2 , . . . ,x n ) (7) 

where i — 1, 2, . . . , n. An alternative to viewing the transformation as a mapping between 
two different spaces is to view it as being correlating the same point in the same space but 
observed from two different coordinate frames of reference which are subject to a similar 
transformation. 

• Coordinate transformations are described as “proper” when they preserve the handed- 
ness (right- or left-handed) of the coordinate system and “improper” when they reverse 
the handedness. Improper transformations involve an odd number of coordinate axes 
inversions through the origin. 

• Inversion of axes may be called improper rotation while ordinary rotation is described 
as proper rotation. 

• Transformations can be active, when they change the state of the observed object (e.g. 
translating the object in space), or passive when they are based on keeping the state of the 
object and changing the state of the coordinate system from which the object is observed. 
Such distinction is based on an implicit assumption of a more general frame of reference 
in the background. 

• Finally, tensor calculus is riddled with conflicting conventions and terminology. In this 
text we will try to use what we believe to be the most common, clear or useful of all of 


these. 
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2 Preliminaries 

2.1 Introduction 

• A tensor is an array of mathematical objects (usually numbers or functions) which 
transforms according to certain rules under coordinates change. In a d-dimensional space, 
a tensor of rarik-n has d n components which may be specified with reference to a given 
coordinate system. Accordingly, a scalar, such as temperature, is a rank-0 tensor with 
(assuming 3D space) 3° = 1 component, a vector, such as force, is a rank-1 tensor with 
3 1 = 3 components, and stress is a rank-2 tensor with 3 2 = 9 components. 

• The term “tensor” was originally derived from the Latin word “tensus” which means 
tension or stress since one of the first uses of tensors was related to the mathematical 
description of mechanical stress. 

• The d n components of a tensor are identified by n distinct integer indices (e.g. i,j,k) 
which are attached, according to the commonly-employed tensor notation, as superscripts 
or subscripts or a mix of these to the right side of the symbol utilized to label the tensor 
(e.g. Aijk, A l i k and A 2fc ). Each tensor index takes all the values over a predefined range 
of dimensions such as 1 to d in the above example of a d-dimensional space. In general, 
all tensor indices have the same range, i.e. they are uniformly dimensioned. 

• When the range of tensor indices is not stated explicitly, it is usually assumed to have 
the values (1,2,3). However, the range must be stated explicitly or implicitly to avoid 
ambiguity. 

• The characteristic property of tensors is that they satisfy the principle of invariance un- 
der certain coordinate transformations. Therefore, formulating the fundamental physical 
laws in a tensor form ensures that they are form-invariant; hence they are objectively- 
representing the physical reality and do not depend on the observer. Having the same 
form in different coordinate systems may be labeled as being “covariant” but this word is 
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also used for a different meaning in tensor calculus as explained in § 2.5.1. 

• “Tensor term” is a product of tensors including scalars and vectors. 

• “Tensor expression” is an algebraic sum (or more generally a linear combination) of 
tensor terms which may be a trivial sum in the case of a single term. 

• “Tensor equality” (symbolized by ‘ = ’) is an equality of two tensor terms and/or expres- 
sions. A special case of this is tensor identity which is an equality of general validity (the 
symbol '=’ may be used for identity as well as for definition). 

• The order of a tensor is identified by the number of its indices (e.g. A l - k is a tensor of 
order 3) which normally identifies the tensor rank as well. However, when contraction (see 
§ 2.6.4) takes place once or more, the order of the tensor is not affected but its rank is 
reduced by two for each contraction operation. 1 

• “Zero tensor” is a tensor whose all components are zero. 

• “Unit tensor” or “unity tensor”, which is usually defined for rank-2 tensors, is a tensor 
whose all elements are zero except the ones with identical values of all indices which are 
assigned the value 1. 

• While tensors of rank-0 are generally represented in a common form of light face non- 

indexed symbols, tensors of rank > 1 are represented in several forms and notations, 

the main ones are the index-free notation, which may also be called direct or symbolic or 

Gibbs notation, and the indicial notation which is also called index or component or tensor 

notation. The first is a geometrically oriented notation with no reference to a particular 

reference frame and hence it is intrinsically invariant to the choice of coordinate systems, 

whereas the second takes an algebraic form based on components identified by indices 

and hence the notation is suggestive of an underlying coordinate system, although being 

a tensor makes it form- invariant under certain coordinate transformations and therefore 

Un the literature of tensor calculus, rank and order of tensors are generally used interchangeably; 
however some authors differentiate between the two as they assign order to the total number of indices, 
including repetitive indices, while they keep rank to the number of free indices. We think the latter is 
better and hence we follow this convention in the present text. 
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it possesses certain invariant properties. The index-free notation is usually identified by 
using bold face symbols, like a and B, while the indicial notation is identified by using 
light face indexed symbols such as a* and B i: . 


2.2 General Rules 

• An index that occurs once in a tensor term is a “free index” . 

• An index that occurs twice in a tensor term is a “dummy” or “bound” index. 

• No index is allowed to occur more than twice in a legitimate tensor term. 2 

• A free index should be understood to vary over its range (e.g. 1 ,n) and hence it 
can be interpreted as saying “for all components represented by the index” . Therefore a 
free index represents a number of terms or expressions or equalities equal to the number 
of allowed values of its range. For example, when i and j can vary over the range 1, . . . , n 
the following expression 

A-i + Bi ( 8 ) 

represents n separate expressions while the following equation 

= B{ (9) 


represents n x n separate equations. 

• According to the “summation convention”, which is widely used in the literature of 
tensor calculus including in the present notes, dummy indices imply summation over their 


2 We adopt this assertion, which is common in the literature of tensor calculus, as we think it is suitable 
for this level. However, there are many instances in the literature of tensor calculus where indices are 
repeated more than twice in a single term. The bottom line is that as long as the tensor expression makes 
sense and the intention is clear, such repetitions should be allowed with no need in our view to take special 
precaution like using parentheses. In particular, the summation convention will not apply automatically 
in such cases although summation on such indices can be carried out explicitly, by using the summation 
symbol or by special declaration of such intention similar to the summation convention. Anyway, in 
the present text we will not use indices repeated more than twice in a single term. 
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range, e.g. for an nD space 


n 


I"/;, = ^17(. - i 1 /: - l 2 /;, + ... + A"n„ 


( 10 ) 


i = 1 


n n 


.vv 


(ii) 


*=1 j=l 


n n n 


e ijk A ij B k = EEE 


( 12 ) 


4=1 j= 1 fc = l 


• When dummy indices do not imply summation, the situation must be clarified by en- 
closing such indices in parentheses or by underscoring or by using upper case letters (with 
declaration of these conventions) or by adding a clarifying comment like “no summation 
on repeated indices” . 

• Tensors with subscript indices, like A i3 , are called covariant, while tensors with super- 
script indices, like A k , are called contravariant. Tensors with both types of indices, like 
A 1 ™ 1 , are called mixed type. More details about this will follow in § 2.5.1. 

• Subscript indices, rather than subscripted tensors, are also dubbed “covariant” and 
superscript indices are dubbed “contravariant”. 

• Each tensor index should conform to one of the variance transformation rules as given 
by Eqs. 20 and 21, i.e. it is either covariant or contravariant. 

• For orthonormal Cartesian coordinate systems, the two variance types (i.e. covariant 
and contravariant) do not differ because the metric tensor is given by the Kronecker delta 
(refer to § 5 and 3.1) and hence any index can be upper or lower although it is common 
to use lower indices in such cases. 

• For tensor invariance, a pair of dummy indices should in general be complementary 
in their variance type, i.e. one covariant and the other contravariant. However, for or- 
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thonormal Cartesian systems the two are the same and hence when both dummy indices 
are covariant or both are contravariant it should be understood as an indication that the 
underlying coordinate system is orthonormal Cartesian if the possibility of an error is 
excluded. 

• As indicated earlier, tensor order is equal to the number of its indices while tensor rank is 
equal to the number of its free indices; hence vectors (terms, expressions and equalities) are 
represented by a single free index and rank- 2 tensors are represented by two free indices. 
The dimension of a tensor is determined by the range taken by its indices. 

• The rank of all terms in legitimate tensor expressions and equalities must be the same. 

• Each term in valid tensor expressions and equalities must have the same set of free 
indices (e.g. i,j, k ). 

• A free index should keep its variance type in every term in valid tensor expressions and 
equations, i.e. it must be covariant in all terms or contravariant in all terms. 

• While free indices should be named uniformly in all terms of tensor expressions and 
equalities, dummy indices can be named in each term independently, e.g. 

4k + B’ t + C'Z (13) 

• A free index in an expression or equality can be renamed uniformly using a different 
symbol, as long as this symbol is not already in use, assuming that both symbols vary 
over the same range, i.e. have the same dimension. 

• Examples of legitimate tensor terms, expressions and equalities: 

4}’ < + C„ = A it - B tJ , a=B j 


(14) 
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• Examples of illegitimate tensor terms, expressions and equalities: 

Bf, Ai + Bij , A* + B j , Ai-B \ A\ = B U (15) 

• Indexing is generally distributive over the terms of tensor expressions and equalities, e.g. 

[A + B], = [A], + [B], (16) 

and 

[A — B] j; <=c [A], = [B], (17) 

• Unlike scalars and tensor components, which are essentially scalars in a generic sense, 
operators cannot in general be freely reordered in tensor terms, therefore 

fh = hf & A r B l = B i A i (18) 

but 

diAi ± AA (19) 

• Almost all the identities in the present notes which are given in a covariant or a con- 
travariant or a mixed form are similarly valid for the other forms unless it is stated other- 
wise. The objective of reporting in only one form is conciseness and to avoid unnecessary 
repetition. 

2.3 Examples of Tensors of Different Ranks 

• Examples of rank-0 tensors (scalars) are energy, mass, temperature, volume and density. 
These are totally identified by a single number regardless of any coordinate system and 
hence they are invariant under coordinate transformations. 
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• Examples of rank-1 tensors (vectors) are displacement, force, electric field, velocity and 
acceleration. These need for their complete identification a number, representing their 
magnitude, and a direction representing their geometric orientation within their space. 
Alternatively, they can be uniquely identified by a set of numbers, equal to the number 
of dimensions of the underlying space, in reference to a particular coordinate system and 
hence this identification is system-dependent although they still have system- invariant 
properties such as length. 

• Examples of rank-2 tensors are Kronecker delta (see § 3.1), stress, strain, rate of strain 
and inertia tensors. These require for their full identification a set of numbers each of 
which is associated with two directions. 

• Examples of rank-3 tensors are the Levi-Civita tensor (see § 3.2) and the tensor of 
piezoelectric moduli. 

• Examples of rank-4 tensors are the elasticity or stiffness tensor, the compliance tensor 
and the fourth-order moment of inertia tensor. 

• Tensors of high ranks are relatively rare in science and engineering. 

• Although rank-0 and rank-1 tensors are, respectively, scalars and vectors, not all scalars 
and vectors (in their generic sense) are tensors of these ranks. Similarly, rank-2 tensors 
are normally represented by matrices but not all matrices represent tensors. 

2.4 Applications of Tensors 

• Tensor calculus is very powerful mathematical tool. Tensor notation and techniques 
are used in many branches of science and engineering such as fluid mechanics, contin- 
uum mechanics, general relativity and structural engineering. Tensor calculus is used for 
elegant and compact formulation and presentation of equations and identities in mathe- 
matics, science and engineering. It is also used for algebraic manipulation of mathematical 
expressions and proving identities in a neat and succinct way (refer to § 4.6). 
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• As indicated earlier, the invariance of tensor forms serves a theoretically and practically 
important role by allowing the formulation of physical laws in coordinate- free forms. 

2.5 Types of Tensors 

• In the following subsections we introduce a number of tensor types and categories and 
highlight their main characteristics and differences. These types and categories are not 
mutually exclusive and hence they overlap in general; moreover they may not be exhaustive 
in their classes as some tensors may not instantiate any one of a complementary set of 
types such as being symmetric or anti-symmetric. 

2.5.1 Covariant and Contravariant Tensors 

• These are the main types of tensor with regard to the rules of their transformation 
between different coordinate systems. 

• Covariant tensors are notated with subscript indices (e.g. Ai) while contravariant tensors 
are notated with superscript indices (e.g. A y ). 

• A covariant tensor is transformed according to the following rule 

A = W A, (20) 

while a contravariant tensor is transformed according to the following rule 

A 1 = ^-A> (21) 

dxi y ’ 

where the barred and unbarred symbols represent the same mathematical object (tensor 
or coordinate) in the transformed and original coordinate systems respectively. 

• An example of covariant tensors is the gradient of a scalar held. 

• An example of contravariant tensors is the displacement vector. 
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• Some tensors have mixed variance type, i.e. they are covariant in some indices and 
contravariant in others. In this case the covariant variables are indexed with subscripts 
while the contravariant variables are indexed with superscripts, e.g. A 1 which is covariant 
in % and contravariant in j. 

• A mixed type tensor transforms covariantly in its covariant indices and contravariantly 
in its contravariant indices, e.g. 


A 1 


n 

m 


dx l dxi dx n 
dx i dx m dx k 



( 22 ) 


• To clarify the pattern of mathematical transformation of tensors, we explain step-by- 
step the practical rules to follow in writing tensor transformation equations between two 
coordinate systems, unbarred and barred, where for clarity we color the symbols of the 
tensor and the coordinates belonging to the unbarred system with blue while we use 
red to mark the symbols belonging to the barred system. Since there are three types 
of tensors: covariant, contravariant and mixed, we use three equations in each step. In 
this demonstration we use rank-4 tensors as examples since this is sufficiently general 
and hence adequate to elucidate the rules for tensors of any rank. The demonstration 
is based on the assumption that the transformation is taking place from the unbarred 
system to the barred system; the same rules should apply for the opposite transformation 
from the barred system to the unbarred system. We use the sign ‘A’ for the equality in 
the transitional steps to indicate that the equalities are under construction and are not 
complete. 

We start by the very generic equations between the barred tensor A and the unbarred 
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tensor A for the three types: 

A = A (covariant) 

A = A (contravariant) (23) 

A = A (mixed) 

We assume that the barred tensor and its coordinates are indexed with ijkl and the 

unbarred are indexed with npqr , so we add these indices in their presumed order and 
position (lower or upper) paying particular attention to the order in the mixed type: 

Aijki = 

A ij kl = 

Since the barred and unbarred tensors are of the same type, as they represent the same 
tensor in two coordinate systems, 3 the indices on the two sides of the equalities should 
match in their position and order. We then insert a number of partial differential operators 
on the right hand side of the equations equal to the rank of these tensors, which is 4 in our 
example. These operators represent the transformation rules for each pair of corresponding 
coordinates one from the barred and one from the unbarred: 

A — 2. Q 2. H 
riijki — d a a a 

Jijkl ° d d d 8 
71 ~ d d d d 

jij — a a a a 
71 ki ~ a a a a 

Now we insert the coordinates of the barred system into the partial differential operators 


-tinpqr 

A n Pqr (25) 

A Tip 

A qr 


npqr 

A np qr ( 24 ) 

Anp 
A qr 


3 Similar basis vectors are assumed. 
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noting that (i) the positions of any index on the two sides should match, i.e. both upper 
or both lower, since they are free indices in different terms of tensor equalities, (ii) a 
superscript index in the denominator of a partial derivative is in lieu of a covariant index 
in the numerator 1 , and (iii) the order of the coordinates should match the order of the 
indices in the tensor: 


a.. -s_ a a a a 

ijkl g x i g x j g x k g x l 

Jijkl dx 1 dxi dx k 8x l 

^ ~ a a a a 

Atj _i_ dx * dxl g g 

kl d d dx k dx 1 


A 

npqr 

j^npqr 


A np 


qr 


( 26 ) 


For consistency, these coordinates should be barred as they belong to the barred tensor; 
hence we add bars: 


-A-ijkl 

a a 

ax'- dxi 

a 

dx k 

a 

dx 1 

A 

npqr 


fiijkl 

° dx 1 dxi 

~ a a 

dx k 

a 

dx 1 

a 

A npqr 

( 27 ) 

A ij 

A kl 

° dx 1 dxi 

~ a a 

a 

dx k 

a 

dx 1 

A Tip 

A qr 



Finally, we insert the coordinates of the unbarred system into the partial differential 
operators noting that (i) the positions of the repeated indices on the same side should 
be opposite, i.e. one upper and one lower, since they are dummy indices and hence the 
position of the index of the unbarred coordinate should be opposite to its position in the 
unbarred tensor, (ii) an upper index in the denominator is in lieu of a lower index in the 
numerator, and (iii) the order of the coordinates should match the order of the indices in 


4 The use of upper indices in the denominator of partial derivatives, which is common in this type of 
equations, is to indicate the fact that the coordinates and their differentials transform contravariantly. 
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the tensor: 



_ dx n dx p dx q dx r a 

dx 1 dx ! dx k dx 1 npqr 


A ijkl 


dx 1 dx q dx k dx l a npqr 
dx n dxP dx q dx r 71 



dx 1 dx q dx q dx r A n P 
dx n dxP dx k dx 1 71 qr 


( 28 ) 


We also replaced the '=’ sign in the final set of equations with the strict equality sign ‘=‘ 
as the equations now are complete. 

• A tensor of m contravariant indices and n covariant indices may be called type (m, n ) 
tensor, e.g. A ? h is a type (1, 2) tensor. When one or both variance types are absent, zero 
is used to refer to the absent type in this notation, e.g. B lk is a type (2, 0) tensor. 

• The covariant and contravariant types of a tensor are linked through the metric tensor 
(refer to § 5). 

• For orthonormal Cartesian systems there is no difference between covariant and con- 
travariant tensors, and hence the indices can be upper or lower. 

• The vectors providing the basis set (not necessarily of unit length or mutually orthogonal) 
for a coordinate system are of covariant type when they are tangent to the coordinate axes, 
and they are of contravariant type when they are perpendicular to the local surfaces of 
constant coordinates. These two sets are identical for orthonormal Cartesian systems. 

• Formally, the covariant and contravariant basis vectors are given respectively by: 

dr 

E, = — & E‘ = Vu* (29) 

du l v ' 

where r is the position vector in Cartesian coordinates and u l is a generalized curvilinear 
coordinate. As indicated already, a superscript in the denominator of partial derivatives 
is equivalent to a subscript in the numerator. 

• In general, the covariant and contravariant basis vectors are not mutually orthogonal 
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or of unit length; however the two sets are reciprocal systems and hence they satisfy the 
following reciprocity relation: 

Ej • E j = Sj (30) 

where Si is the Kronecker delta (refer to § 3.1). 

• A vector can be represented either by covariant components with contravariant coordi- 
nate basis vectors or by contravariant components with covariant coordinate basis vectors. 
For example, a vector A can be expressed as 

A = A, E' or A = A% (31) 

where E' and E; are the contravariant and covariant basis vectors respectively. The use of 
the covariant or contravariant form of the vector representation is a matter of choice and 
convenience. 

• More generally, a tensor of any rank (> 1) can be represented covariantly using con- 
travariant basis tensors of that rank, or contravariantly using covariant basis tensors, or 
in a mixed form using a mixed basis of opposite type. For example, a rank-2 tensor A can 
be written as: 

A = ,1,;E E ; = A ij EiEj = A/E% (32) 

where E'E J . E,-E y and E'E y are dyadic products (refer to § 2.6.3). 

2.5.2 True and Pseudo Tensors 

• These are also called polar and axial tensors respectively although it is more common 
to use the latter terms for vectors. Pseudo tensors may also be called tensor densities. 

• True tensors are proper (or ordinary) tensors and hence they are invariant under co- 
ordinate transformations, while pseudo tensors are not proper tensors since they do not 
transform invariantly as they acquire a minus sign under improper orthogonal transfer- 
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mations which involve inversion of coordinate axes through the origin with a change of 
system handedness. 

• Because true and pseudo tensors have different mathematical properties and represent 
different types of physical entities, the terms of consistent tensor expressions and equations 
should be uniform in their true and pseudo type, i.e. all terms are true or all are pseudo. 

• The direct product (refer to § 2.6.3) of even number of pseudo tensors is a proper tensor, 
while the direct product of odd number of pseudo tensors is a pseudo tensor. The direct 
product of true tensors is obviously a true tensor. 

• The direct product of a mix of true and pseudo tensors is a true or pseudo tensor 
depending on the number of pseudo tensors involved in the product as being even or odd 
respectively. 

• Similar rules to those of direct product apply to cross products (including curl operations) 
involving tensors (usually of rank-1) with the addition of a pseudo factor for each cross 
product operation. This factor is contributed by the permutation tensor (see § 3.2) which 
is implicit in the definition of the cross product (see Eqs. 121 and 146). 

• In summary, what determines the tensor type (true or pseudo) of the tensor terms in- 
volving direct 5 and cross products is the parity of the multiplicative factors of pseudo type 
plus the number of cross product operations involved since each cross product contributes 
an e factor. 

• Examples of true scalars are temperature, mass and the dot product of two polar or two 
axial vectors, while examples of pseudo scalars are the dot product of an axial vector and 
a polar vector and the scalar triple product of polar vectors. 

• Examples of polar vectors are displacement and acceleration, while examples of axial 

vectors are angular velocity and cross product of polar vectors in general (including curl 

operation on polar vectors) due to the involvement of the permutation symbol which is 

5 Inner product (see § 2.6.5) is the result of a direct product operation followed by a contraction (see 
§ 2.6.4) and hence it is a direct product in this context. 
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a pseudo tensor (refer to § 3.2). The essence of this distinction is that the direction of a 
pseudo vector depends on the observer choice of the handedness of the coordinate system 
whereas the direction of a proper vector is independent of such choice. 

• Examples of proper tensors of rank-2 are stress and rate of strain tensors, while examples 
of pseudo tensors of rank-2 are direct products of two vectors: one polar and one axial. 

• Examples of proper tensors of higher ranks are piezoelectric moduli tensor (rank-3) 
and elasticity tensor (rank-4), while examples of pseudo tensors of higher ranks are the 
permutation tensor of these ranks. 


2.5.3 Absolute and Relative Tensors 

• Considering an arbitrary transformation from a general coordinate system to another, a 
relative tensor of weight w is defined by the following tensor transformation: 


Aij...k 


dx w dx 1 dx j dx k dx d dx e dx^ b _ c 

dx dx a dx b dx c dx 1 dx m dx n d 


(33) 


where |||| is the Jacobian of the transformation between the two systems. When w = 0 
the tensor is described as an absolute or true tensor, while when w — — 1 the tensor is 
described as a pseudo tensor. When w — 1 the tensor may be described as a tensor 
density. 6 

• As indicated earlier, a tensor of m contravariant indices and n covariant indices may be 
called type (m,n). This may be generalized to include the weight as a third entry and 
hence the type of the tensor is identified by (m,n,w). 

• Relative tensors can be added and subtracted if they are of the same variance type and 
have the same weight; the result is a tensor of the same type and weight. Also, relative 
tensors can be equated if they are of the same type and weight. 

• Multiplication of relative tensors produces a relative tensor whose weight is the sum of 


6 Some of these labels are used differently by different authors. 
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the weights of the original tensors. Hence, if the weights are added up to a non- zero value 
the result is a relative tensor of that weight; otherwise it is an absolute tensor. 

2.5.4 Isotropic and Anisotropic Tensors 

• Isotropic tensors are characterized by the property that the values of their components 
are invariant under coordinate transformation by proper rotation of axes. In contrast, the 
values of the components of anisotropic tensors are dependent on the orientation of the 
coordinate axes. Notable examples of isotropic tensors are scalars (rank-0), the vector 0 
(rank-1), Kronecker delta (rank-2) and Levi-Civita tensor e l3 k (rank-3). Many tensors 
describing physical properties of materials, such as stress and magnetic susceptibility, are 
anisotropic. 

• Direct and inner products of isotropic tensors are isotropic tensors. 

• The zero tensor of any rank is isotropic; therefore if the components of a tensor vanish 
in a particular coordinate system they will vanish in all properly and improperly rotated 
coordinate systems.' Consequently, if the components of two tensors are identical in a 
particular coordinate system they are identical in all transformed coordinate systems. 

• As indicated, all rank-0 tensors (scalars) are isotropic. Also, the zero vector, 0, of any 
dimension is isotropic; in fact it is the only rank-1 isotropic tensor. 

2.5.5 Symmetric and Anti-symmetric Tensors 

• These types of tensor apply to high ranks only (rank > 2). Moreover, these types are 
not exhaustive, even for tensors of ranks > 2, as there are high-rank tensors which are 
neither symmetric nor anti-symmetric. 


7 For improper rotation, this is more general than being isotropic. 
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• A rank- 2 tensor A l3 is symmetric iff for all i and j 


A .. — A.. 

-^Jl Sl-lJ 


(34) 


and anti-symmetric or skew-symmetric iff 


Aji — 


-Aij 


(35) 


Similar conditions apply to contravariant type tensors (refer also to the following). 
• A rank -n tensor is symmetric in its two indices ij and i[ iff 


A. . . . — A. . . 

■ rv l\...ll...lj...ln - rL ll...lj...ll...ln 


and anti-symmetric or skew-symmetric in its two indices ij and %i iff 


(36) 


A. ... — _A. . . . 


(37) 


• Any rank- 2 tensor A l3 can be synthesized from (or decomposed into) a symmetric part 
(marked with round brackets enclosing the indices) and an anti-symmetric part Auj] 
(marked with square brackets) where 


Aij A(jj) + A[jj], 


A(ij) 2 (Aij + Aji ) & 


A[ij] 2 ( A%j Aji) (38) 


• A rank-3 tensor A ljk can be symmetrized by 

A(ijk) yy (A ijk A ki j Aj k i -\- Ai k j -\- Aji k -\- A k ji) 


(39) 
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and anti-symmetrized by 


A\jjk\ ( A 3 k T A-kij + Ajki 


A-ikj Ajik Akji) 


(40) 


• A rank-n tensor can be symmetrized by 


A 


(il-.in) 


— (sum of all even & odd permutations of indices i’s) 
n! 


(41) 


and anti-symmetrized by 

Au i i 1 = — r (sum of all even permutations minus sum of all odd permutations) (42) 
n\ 

• For a symmetric tensor A tJ and an anti-symmetric tensor B l i (or the other way around) 
we have 

AijB ij = 0 (43) 

• The indices whose exchange defines the symmetry and anti-symmetry relations should 
be of the same variance type, i.e. both upper or both lower. 

• The symmetry and anti-symmetry characteristic of a tensor is invariant under coordinate 
transformation. 

• A tensor of high rank (> 2) may be symmetrized or anti-symmetrized with respect to 
only some of its indices instead of all of its indices, e.g. 

A(ijy. 2 T Ajik') A — (Aijk Ajik ) (44) 


• A tensor is totally symmetric iff 


A. 


...lr 


= A 


0 il-.in ) 


(45) 



2.6 Tensor Operations 


and totally anti-symmetric iff 


A. ■ — A r- • 

■ ri -ll...ln -^*-[11... In 


(46) 


• For a totally skew-symmetric tensor (i.e. anti-symmetric in all of its indices), nonzero 
entries can occur only when all the indices are different. 


2.6 Tensor Operations 

• There are many operations that can be performed on tensors to produce other tensors 
in general. Some examples of these operations are addition/subtraction, multiplication 
by a scalar (rank-0 tensor), multiplication of tensors (each of rank > 0), contraction and 
permutation. Some of these operations, such as addition and multiplication, involve more 
than one tensor while others are performed on a single tensor, such as contraction and 
permutation. 

• In tensor algebra, division is allowed only for scalars, hence if the components of an 
indexed tensor should appear in a denominator, the tensor should be redefined to avoid 
this, e.g. Bi = 

2.6.1 Addition and Subtraction 

• Tensors of the same rank and type (covariant /contravariant /mixed and true/pseudo) 
can be added algebraically to produce a tensor of the same rank and type, e.g. 

a = b + c (47) 


Ai — Bi — Ci 


(48) 
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A) = B) + C) (49) 

• The added/subtracted terms should have the same indicia! structure with regard to 
their free indices, as explained in § 2.2, hence A’- k and B J ik cannot be added or subtracted 
although they are of the same rank and type, but and B' l - k can be added and sub- 
tracted. 

• Addition of tensors is associative and commutative: 

(A + B) + C = A + (B + C) (50) 

A + B = B + A (51) 


2.6.2 Multiplication by Scalar 

• A tensor can be multiplied by a scalar, which generally should not be zero, to produce 
a tensor of the same variance type and rank, e.g. 

A = (52) 


where a is a non- zero scalar. 

• As indicated above, multiplying a tensor by a scalar means multiplying each component 
of the tensor by that scalar. 

• Multiplication by a scalar is commutative, and associative when more than two factors 


are involved. 
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2.6.3 Tensor Multiplication 

• This may also be called outer or exterior or direct or dyadic multiplication, although 
some of these names may be reserved for operations on vectors. 

• On multiplying each component of a tensor of rank r by each component of a tensor of 
rank k, both of dimension m, a tensor of rank (r + k) with m r+k components is obtained 
where the variance type of each index (covariant or contravariant) is preserved, e.g. 

AiBj = Cn (53) 

A ,s B a = C*, (54) 

• The outer product of a tensor of type (m, n) by a tensor of type (p, q ) results in a tensor 
of type (m + p,n + q). 

• Direct multiplication of tensors may be marked by the symbol ®, mostly when using 
symbolic notation for tensors, e.g. A ® B. However, in the present text no symbol will be 
used for the operation of direct multiplication. 

• Direct multiplication of tensors is not commutative. 

• The outer product operation is distributive with respect to the algebraic sum of tensors: 

A (B ± C) = AB ± AC & (B ± C) A = BA ± CA (55) 

• Multiplication of a tensor by a scalar (refer to § 2.6.2) may be regarded as a special case 
of direct multiplication. 

• The rank-2 tensor constructed as a result of the direct multiplication of two vectors is 
commonly called dyad. 

• Tensors may be expressed as an outer product of vectors where the rank of the resultant 
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product is equal to the number of the vectors involved (e.g. 2 for dyads and 3 for triads). 

• Not every tensor can be synthesized as a product of lower rank tensors. 

• In the outer product, it is understood that all the indices of the involved tensors have 
the same range. 

• The outer product of tensors yields a tensor. 

2.6.4 Contraction 

• Contraction of a tensor of rank > 1 is to make two free indices identical, by unifying 
their symbols, and perform summation over these repeated indices, e.g. 


• Contraction results in a reduction of the rank by 2 since it implies the annihilation of 
two free indices. Therefore, the contraction of a rank-2 tensor is a scalar, the contraction 
of a rank-3 tensor is a vector, the contraction of a rank-4 tensor is a rank-2 tensor, and so 
on. 

• For general non-Cartesian coordinate systems, the pair of contracted indices should be 
different in their variance type, i.e. one upper and one lower. Hence, contraction of a 
mixed tensor of type (m,n) will, in general, produce a tensor of type (m — l,n — 1). 

• A tensor of type (p, q ) can have p x q possible contractions, i.e. one contraction for each 
pair of lower and upper indices. 

• A common example of contraction is the dot product operation on vectors which can be 
regarded as a direct multiplication (refer to § 2.6.3) of the two vectors, which results in a 
rank-2 tensor, followed by a contraction. 


Aj contraction A\ 


(56) 


.4)) contraction on jl A 


(57) 
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• In matrix algebra, taking the trace (summing the diagonal elements) can also be consid- 
ered as contraction of the matrix, which under certain conditions can represent a rank-2 
tensor, and hence it yields the trace which is a scalar. 

• Applying the index contraction operation on a tensor results into a tensor. 

• Application of contraction of indices operation on a relative tensor (see § 2.5.3) produces 
a relative tensor of the same weight as the original tensor. 

2.6.5 Inner Product 

• On taking the outer product (refer to § 2.6.3) of two tensors of rank > 1 followed by a 
contraction on two indices of the product, an inner product of the two tensors is formed. 
Hence if one of the original tensors is of rank-m and the other is of rank-n, the inner 
product will be of rank-(m + n — 2). 

• The inner product operation is usually symbolized by a single dot between the two 
tensors, e.g. A • B, to indicate contraction following outer multiplication. 

• In general, the inner product is not commutative. When one or both of the tensors 
involved in the inner product are of rank > 1 the order of the multiplicands does matter. 

• The inner product operation is distributive with respect to the algebraic sum of tensors: 

A-(B±C) = A- B±A-C & (B±C)-A = B- A±C-A (58) 

• As indicated before (see § 2.6.4), the dot product of two vectors is an example of the 
inner product of tensors, i.e. it is an inner product of two rank-1 tensors to produce a 
rank-0 tensor: 

[ab]/ = dilP contraction^ a • b = (59) 

• Another common example (from linear algebra) of inner product is the multiplication of 
a matrix (representing a rank-2 tensor assuming certain conditions) by a vector (rank-1 
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tensor) to produce a vector, e.g. 


[Ab]^ = Aijb k contraction on jk [A • b]. = A t] P 


( 60 ) 


The multiplication of two n x n matrices is another example of inner product (see Eq. 
119). 

• For tensors whose outer product produces a tensor of rank > 2, various contraction 
operations between different sets of indices can occur and hence more than one inner 
product, which are different in general, can be defined. Moreover, when the outer product 
produces a tensor of rank > 3 more than one contraction can take place simultaneously. 

• There are more specialized types of inner product; some of which may be defined dif- 
ferently by different authors. For example, a double inner product of two rank-2 tensors, 
A and B, may be defined and denoted by double vertically- or horizontally- aligned dots 
(e.g. A: B or A • -B) to indicate double contraction taking place between different pairs 
of indices. An instance of these types is the inner product with double contraction of two 
dyads which is commonly defined by 8 


and the result is a scalar. The single dots in the right hand side of the last equation 
symbolize the conventional dot product of two vectors. Some authors may define a different 
type of double-contraction inner product of two dyads, symbolized by two horizontally- 
aligned dots, which may be called a “transposed contraction” , and is given by 


ab : cd = (a • c) (b • d) 


( 61 ) 


ab ■ cd = ab : dc = (a ■ d) (b ■ c) 


( 62 ) 


8 It is also defined differently by some authors. 
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where the result is also a scalar. However, different authors may have different conventions 
and hence one should be vigilant about such differences. 

• For two rank- 2 tensors, the aforementioned double-contraction inner products are simi- 
larly defined as in the case of two dyads: 

A: B = AijBij & A • B = AijBji (63) 

• Inner products with higher multiplicity of contractions are similarly defined, and hence 
can be regarded as trivial extensions of the inner products with lower contraction multi- 
plicities. 

• The inner product of tensors produces a tensor because the inner product is an outer 
product operation followed by an index contraction operation and both of these operations 
on tensors produce tensors. 

2.6.6 Permutation 

• A tensor may be obtained by exchanging the indices of another tensor, e.g. transposition 
of rank-2 tensors. 

• Tensor permutation applies only to tensors of rank > 2. 

• The collection of tensors obtained by permuting the indices of a basic tensor may be 
called isomers. 

2.7 Tensor Test: Quotient Rule 

• Sometimes a tensor-like object may be suspected for being a tensor; in such cases a test 
based on the “quotient rule” can be used to clarify the situation. According to this rule, if 
the inner product of a suspected tensor with a known tensor is a tensor then the suspect 
is a tensor. In more formal terms, if it is not known if A is a tensor but it is known that 



2.7 Tensor Test: Quotient Rule 


35 


B and C are tensors; moreover it is known that the following relation holds true in all 
rotated (properly-transformed) coordinate frames: 

Apq...k...mBij ...k...n C pq...mij...n (64) 

then A is a tensor. Here, A, B and C are respectively of ranks m, n and (m + n — 2), due 
to the contraction on k which can be any index of A and B independently. 

• Testing for being a tensor can also be done by applying first principles through direct 
substitution in the transformation equations. However, using the quotient rule is generally 
more convenient and requires less work. 

• The quotient rule may be considered as a replacement for the division operation which 
is not defined for tensors. 
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3 5 and e Tensors 

• These tensors are of particular importance in tensor calculus due to their distinctive 
properties and unique transformation attributes. They are numerical tensors with fixed 
components in all coordinate systems. The first is called Kronecker delta or unit ten- 
sor, while the second is called Levi-Civita 9 , permutation, anti-symmetric and alternating 
tensor. 

• The 6 and e tensors are conserved under coordinate transformations and hence they are 
the same for all systems of coordinate. 10 


3.1 Kronecker 5 

• This is a rank-2 symmetric tensor in all dimensions, i.e. 

Sij = Sji (i,j = 1,2,. . . ,ra) (65) 


Similar identities apply to the contravariant and mixed types of this tensor. 

• It is invariant in all coordinate systems, and hence it is an isotropic tensor. * 11 

• It is defined as: 


Jij 


1 

0 


(* = j) 
(* + j) 


( 66 ) 


9 This name is usually used for the rank-3 tensor. Also some authors distinguish between the permuta- 
tion tensor and Levi-Civita tensor even for rank-3. Moreover, some of the common labels and descriptions 
of e are more specific to rank-3. 

10 For the permutation tensor, the statement applies to proper coordinate transformations. 

11 In fact it is more general than isotropic as it is invariant even under improper coordinate transfor- 
mations. 
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and hence it can be considered as the identity matrix, e.g. for 3D 


1 

on 

812 

On 

CO 

1 


1 

0 

1 

0 

5 2 i 

822 

823 

= 

0 

1 

0 

CO 

832 

833 


0 

0 

1 


(67) 


• Covariant, contravariant and mixed type of this tensor are the same, that is 


8) = 8? = 5 ij = Sij 


( 68 ) 


3.2 Permutation e 

• This is an isotropic tensor. It has a rank equal to the number of dimensions; hence, a 
rank -n permutation tensor has n n components. 

• It is totally anti-symmetric in each pair of its indices, i.e. it changes sign on swapping 
any two of its indices, that is 


The reason is that any exchange of two indices requires an even/odd number of single- 
step shifts to the right of the first index plus an odd/even number of single-step shifts to 
the left of the second index, so the total number of shifts is odd and hence it is an odd 
permutation of the original arrangement. 

• It is a pseudo tensor since it acquires a minus sign under improper orthogonal transfor- 
mation of coordinates (inversion of axes with possible superposition of rotation). 

• Definition of rank- 2 e (e^-): 


Cl2 — 1, 


621 — — 1 & £ii — 622 — 0 


( 70 ) 
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• Definition of rank-3 e (e^-fc): 


^ijk 


= < 


1 

-1 

0 


(■ i,j , k is even permutation of 1,2,3) 
(■ i,j , k is odd permutation of 1,2,3) 
(repeated index) 


(71) 


• The definition of rank-n e (e^...^) is similar to the definition of rank-3 e considering 
index repetition and even or odd permutations of its indices (ii,i 2 , • • • , i n ) corresponding 
to (1, 2, • • • , n ) , that is 


Uli 2 —in ~ \ 


1 

-1 

0 


[(*i, i 2 , • • • , i n ) is even permutation of (1,2,..., n)] 

[(*!, i 2 , • • • , i n ) is odd permutation of (1,2,..., n)] (72) 

[repeated index] 


• e may be considered a contravariant relative tensor of weight +1 or a covariant relative 
tensor of weight —1. Hence, in 2D, 3D and nD spaces respectively we have: 


e .. — 

7 C 


= e ijk 




(73) 

(74) 

(75) 


3.3 Useful Identities Involving 8 or/and e 

3.3.1 Identities Involving 5 

• When an index of the Kronecker delta is involved in a contraction operation by repeating 
an index in another tensor in its own term, the effect of this is to replace the shared index 
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in the other tensor by the other index of the Kronecker delta, that is 




( 76 ) 


In such cases the Kronecker delta is described as the substitution or index replacement 
operator. Hence, 


OijOjk if' 


ik 


( 77 ) 


Similarly, 


Oijfijhvki ^ik^ki ^ 


( 78 ) 


where n is the space dimension. 

• Because the coordinates are independent of each other: 


dxj 

dx i 


dyXi = Xij = 5, 


jij 


( 79 ) 


Hence, in an nD space we have 


diXi = 5a = n 


( 80 ) 


• For orthonormal Cartesian systems: 


dx l 

dxi 


dx j 

dx i 


Sij = S ij 


( 81 ) 


• For a set of orthonormal basis vectors in orthonormal Cartesian systems: 




->IJ 


( 82 ) 


• The double inner product of two dyads formed by orthonormal basis vectors of an 
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orthonormal Cartesian system is given by: 




(83) 


3.3.2 Identities Involving e 

• For rank-3 e: 


eijk = ekij = tjki = ~tikj = —ejik = -ekji (sense of cyclic order) (84) 


These equations demonstrate the fact that rank-3 e is totally anti-symmetric in all of its 
indices since a shift of any two indices reverses the sign. This also reflects the fact that 
the above tensor system has only one independent component. 

• For rank- 2 e: 

Cj = (J ~ i) (85) 

• For rank-3 e: 

Cjfc = \(j ~ ( k -*)(*“ j) (86) 

• For rank-4 e: 


tijU = 4 0 - i) ( k - i) (l - i)(k-j)(l - j) (( - k) 


(87) 


For rank-n e: 


c aia2---a n 


n—1 

n 


■\ ( a i a *) 


i= 1 L j=i+l 


5(n - 1) 


1 1 (°b 


( 88 ) 
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where S(n — 1) is the super-factorial function of (n — 1) which is defined as 

k 

S(k) = JJi! = 1! -2! ■ ... • k\ (89) 

i = 1 

A simpler formula for rank -n e can be obtained from the previous one by ignoring the 
magnitude of the multiplication factors and taking only their signs, that is 


Caia2---a„ 


n a (%■ ~ a *) 

l<i<7<n 



(90) 


where 


• For rank-n e: 


( 



+i 

(k > 0) 


a(k) = < 

-l 

(k < 0) 

(91) 


0 

(k = 0) 



"in ^ill2-"in 

(92) 


because this is the sum of the squares of e ili2 ... in over all the permutations of n different 
indices which is equal to n\ where the value of e of each one of these permutations is either 
+1 or —1 and hence in both cases their square is 1. 

• For a symmetric tensor A ]k : 

CijkAjk = 0 (93) 

because an exchange of the two indices of A jk does not affect its value due to the symmetry 
whereas a similar exchange in these indices in e ljk results in a sign change; hence each term 
in the sum has its own negative and therefore the total sum will vanish. 


^ ijkAiAj 6ijkAiA k 6ij k Aj A k 0 


(94) 
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because, due to the commutativity of multiplication, an exchange of the indices in kl’s will 
not affect the value but a similar exchange in the corresponding indices of will cause 
a change in sign; hence each term in the sum has its own negative and therefore the total 
sum will be zero. 

• For a set of orthonormal basis vectors in a 3D space with a right-handed orthonormal 
Cartesian coordinate system: 

X Gj (95) 

X G/dj ^ ijk (96) 

3.3.3 Identities Involving 5 and e 


^ijk0lid2j0sk ~ Cl23 — 1 


• For rank- 2 e: 




&ik &il 
& jk &jl 


&ik&jl ^il^jk 


(97) 


(98) 


tiltkl ~ 


&ik 


(99) 


^-ij € ij 2 


(100) 
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For rank-3 e: 




^ il $ im &in 
djl Vjm Ojn 
$kl &km &kn 


jm^kn^^im,^ j kl~\~3i n 5 jiS km ^il^jn^km ^im^jl^kn ^in^jm^kl 


( 101 ) 


^ijk^lmk 




( 102 ) 


&il &im 

$jl 3jm 

The last identity is very useful in manipulating and simplifying tensor expressions and 
proving vector and tensor identities. 


^ijk^ljk 2 5u 


(103) 


^ij k^ijk 2 Sn 6 


(104) 


since the rank and dimension of e are the same, which is 3 in this case. 
• For rank-n e: 





5 - • • 

■ 5- ■ 

'nil 

U 11J2 

u lljn 

'nii 

^ni2 

^nin 

injl 

Vinj2 

5i j 

L njn 


(105) 


According to Eqs. 71 and 76: 


^ ijk^ij ^-ijk&ik ^ijk^jk 0 


( 106 ) 
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3.4 Generalized Kronecker delta 

• The generalized Kronecker delta is defined by: 

[(ji • • ■ j n ) is even permutation of (A . . . i n )] 

[(j\ . . . j n ) is odd permutation of {i\ . . . i„)\ (107) 

[repeated j’s] 

It can also be defined by the following n x n determinant: 


r*l ...In 


1 
-1 
0 


CZi . ..In 

U jl-3n 



*2 •' 

■■ 5 ? 

Jn 


% •' 

■ 

Jn 


*2 •• 

■ 

Jn 


(108) 


where the <5* entries in the determinant are the normal Kronecker delta as defined by Eq. 

66 . 

• Accordingly, the relation between the rank-n e and the generalized Kronecker delta in 
an nD space is given by: 


el 2...n 

tiii 2 ...i n °i 1 i 2 ...i n 


& 




= <5i : 


l\l 2 ...l n 
2 .. .72 


(109) 


Hence, the permutation tensor e may be considered as a special case of the generalized 
Kronecker delta. Consequently the permutation symbol can be written as an n x n deter- 
minant consisting of the normal Kronecker deltas. 

• If we define 


cij cijk 
°lm ~ ° Imk 


( 110 ) 
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then Eq. 102 will take the following form: 



s‘A - s'jj 


(in) 


Other identities involving S and e can also be formulated in terms of the generalized 
Kronecker delta. 

• On comparing Eq. 105 with Eq. 108 we conclude 


__ 11... J. 

U h-in 


'-n-u 


( 112 ) 
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4 Applications of Tensor Notation and Techniques 

4.1 Common Definitions in Tensor Notation 

• The trace of a matrix A representing a rank- 2 tensor is: 

tr (A) = An (113) 





( 118 ) 
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It should be noticed that here we are using matrix notation. The multiplication operation, 
according to the symbolic notation of tensors, should be denoted by a dot between the 
tensor and the vector, i.e. Ab. 1 

• The multiplication of two n x n matrices A and B as defined in linear algebra is: 

[AB] ifc = AijBjk (119) 

Again, here we are using matrix notation; otherwise a dot should be inserted between the 
two matrices. 

• The dot product of two vectors is: 


A • B =5 ij A i B j = AiBi (120) 

The readers are referred to § 2.6.5 for a more general definition of this type of product 
that includes higher rank tensors. 

• The cross product of two vectors is: 


[A x B] i eij k AjB k 


( 121 ) 


The scalar triple product of three vectors is: 


A - (B x C) = 


A\ A 2 A 3 
B\ B 2 B 3 

C i c 2 c 3 


eijkAiBjC k 


( 122 ) 


12 The matrix multiplication in matrix notation is equivalent to a dot product operation in tensor 
notation. 
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• The vector triple product of three vectors is: 

[A x (B x C)] i = 6ijkCkimAjBiC m (123) 

4.2 Scalar Invariants of Tensors 

• In the following we list and write in tensor notation a number of invariants of low rank 
tensors which have special importance due to their widespread applications in vector and 
tensor calculus. All These invariants are scalars. 

• The value of a scalar (rank-0 tensor), which consists of a magnitude and a sign, is 
invariant under coordinate transformation. 

• An invariant of a vector (rank-1 tensor) under coordinate transformations is its magni- 
tude, i.e. length (the direction is also invariant but it is not scalar!). 13 

• The main three independent scalar invariants of a rank-2 tensor A under change of basis 
are: 

/ = tr (A) = An (124) 

II = tr (A 2 ) = AijAji (125) 

III = tr (A 3 ) = AijAjkAki (126) 

• Different forms of the three invariants of a rank-2 tensor A, which are also widely used, 
are: 

h = I = A, (127) 


13 In fact the magnitude alone is invariant under coordinate transformations even for pseudo vectors 
because it is a scalar. 
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h = \ (/ 2 - //) = \ (AuAjj - AijAji) (128) 

I 3 = det (A) = — (/ 3 — 31 II + 2 III) = —^ijktpqrAipAjqAkr (129) 

• The invariants I, II and III can similarly be defined in terms of the invariants Ji, 0 
and J 3 as follow: 

I = h (130) 

II = if - 2 1 2 (131) 

III = if - 3hh + 3/ s (132) 

• Since the determinant of a matrix representing a rank-2 tensor is invariant, then if the 


determinant vanishes in one coordinate system it will vanish in all coordinate systems and 
vice versa. Consequently, if a rank- 2 tensor is invertible in a particular coordinate system, 
it is invertible in all coordinate systems. 

• Ten joint invariants between two rank- 2 tensors, A and B, can be formed; these are: 
tr(A), tr(B), tr(A 2 ), tr(B 2 ), tr(A 3 ), tr(B 3 ), tr(A-B), tr(A 2 -B), tr(A-B 2 ) and 
tr (A 2 -B 2 ). 

4.3 Common Differential Operations in Tensor Notation 

• Here we present the most common differential operations as defined by tensor notation. 
These operations are mostly based on the various types of interaction between the vector 
differential operator nabla V with tensors of different ranks as well as interaction with 
other types of operation like dot and cross products. 
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• V is essentially a spatial partial differential operator defined in Cartesian coordinate 
systems by: 

v - = k < 133 > 

The definition of V in some non-Cartesian systems will be given in § 4.3.2. 


4.3.1 Cartesian System 

• The gradient of a differentiable scalar function of position / is a vector given by: 


[V/], = V ?: / 


d£ 

dxi 


dj = f,i 


(134) 


• The gradient of a differentiable vector function of position A (which is the outer product, 
as defined in § 2.6.3, between the V operator and the vector) is a rank-2 tensor defined 
by: 

[VA] y = diAj (135) 

• The gradient operation is distributive but not commutative or associative: 


V (/ + h) = V/ + Vh 


(136) 


v/v/v 


(137) 


(V/) h ± V (fh) 


( 138 ) 


where / and h are differentiable scalar functions of position. 
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• The divergence of a differentiable vector A is a scalar given by: 


V- A 


. 34 

>V dxj 


dAj 

dxi 


VjA i — diAi — A iti 


(139) 


The divergence operation can also be viewed as taking the gradient of the vector followed 
by a contraction. Hence, the divergence of a vector is invariant because it is the trace of 
a rank-2 tensor. 14 

• The divergence of a differentiable rank-2 tensor A is a vector defined in one of its forms 
by: 

[V ■ A], = djAjt (140) 


and in another form by 

[V ■ A] . = diAji (141) 

These two different forms can be given, respectively, in symbolic notation by: 


V ■ A & V • A T 


(142) 


where A 1 is the transpose of A. More generally, the divergence of a tensor of rank n > 2, 
which is a tensor of rank- (n — 1), can be defined in several forms, which are different in 
general, depending on the combination of the contracted indices. 

• The divergence operation is distributive but not commutative or associative: 


V-(A + B) = V- A + V- B 


(143) 


V-A# A-V 


( 144 ) 


14 


It may also be argued that the divergence of a vector is a scalar and hence it is invariant. 
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V ■ (/A) ^ V/ • A (145) 

where A and B are differentiable tensor functions of position. 

• The curl of a differentiable vector A is a vector given by: 


[V x A]^ 6ijk 


d A k 
dxj 


^ijk-^-k.j 


(146) 


• The curl operation may be generalized to tensors of rank > 1, and hence the curl of a 
differentiable rank- 2 tensor A can be defined as a rank- 2 tensor given by: 


[V X A]y 


r) A 


z 'imn KJ m 1 x nj 


• The curl operation is distributive but not commutative or associative: 


(147) 


Vx(A + B)-VxA + VxB 


(148) 


VxA/Ax V 


(149) 


V x (A x B) / (V x A) x B 


(150) 


• The Laplacian scalar operator, also called the harmonic operator, acting on a differen- 
tiable scalar / is given by: 


A/ 


V 2 / 


Jlj 


d 2 f 

dxidxj 


d 2 f 

dxidxi 


V«/ = d uf = /,, 


(151) 


The Laplacian operator acting on a differentiable vector A is defined for each component 
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of the vector similar to the definition of the Laplacian acting on a scalar, that is 

( 152 ) 

• The following scalar differential operator is commonly used in science (e.g. in fluid 
dynamics): 

A • V = AiVi = A^- = Adi (153) 

OXi 

where A is a vector. As indicated earlier, the order of Ai and <9., should be respected. 

• The following vector differential operator also has common applications in science: 

[A x V] ■ = eijkAjdk (154) 

• The differentiation of a tensor increases its rank by one, by introducing an extra covariant 
index, unless it implies a contraction in which case it reduces the rank by one. Therefore 
the gradient of a scalar is a vector and the gradient of a vector is a rank- 2 tensor (<9jA,), 
while the divergence of a vector is a scalar and the divergence of a rank-2 tensor is a vector 
( djAji or d iAji). This may be justified by the fact that V is a vector operator. On the 
other hand the Laplacian operator does not change the rank since it is a scalar operator; 
hence the Laplacian of a scalar is a scalar and the Laplacian of a vector is a vector. 

4.3.2 Other Coordinate Systems 

• For completeness, we define here some differential operations in the most commonly 
used non-Cartesian coordinate systems, namely cylindrical and spherical systems, as well 
as general orthogonal coordinate systems. 

• We can use indexed generalized coordinates like q \ , q 2 and q% for the cylindrical coor- 
dinates (p, 4>, z) and the spherical coordinates (r, d,0). However, for more clarity at this 
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level and to follow the more conventional practice, we use the coordinates of these systems 
as suffixes in place of the indices used in the tensor notation. 15 

• For the cylindrical system identified by the coordinates (p, </>, z ) with an orthonormal 
basis vectors e p , e ^ and e 2 : lb 
The V operator is: 


V 


Qpdp T &<j> dfp T e z d z 
P 


(155) 


The Laplacian operator is: 


V" — d pp + -dp H — ^ 


P A 


The gradient of a differentiable scalar / is: 


d/ . 10/ W 

1 dp ^ pdcj) z dz 


The divergence of a differentiable vector A is: 

V- A = - 
P 

The curl of a differentiable vector A is: 


d (pAp) dA ^ d (pA z ) 

dp d(f) dz 


V x A 


1 

P 


Qp P^<j> ® z 

a_ d_ d_ 

dp d(j> dz 


Ap pA(f) A z 


(156) 


(157) 


(158) 


(159) 


For plane polar coordinate systems, these operators and operations can be obtained by 
dropping the z components or terms from the cylindrical form of the above operators and 


15 There is another reason that is these are physical components not covariant or contravariant. 

16 It should be obvious that since p, <fi and z are specific coordinates and not variable indices, the 
summation convention does not apply. 
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operations. 

• For the spherical system identified by the coordinates (r, 9, (p) with an orthonormal basis 
vectors e r , eg and e^: 1 ' 

The V operator is: 

V = e r d r + eg-dg + T — -8$ (160) 

r r sm 9 

The Laplacian operator is: 


V 2 


drr H d r ~\ T^dgg + 


cos 9 
r 2 sin 9 


dg + 



(161) 


The gradient of a differentiable scalar / is: 


df, ±df 1 df 

r dr 9 r d9 <b r sin 9 dd 


(162) 


The divergence of a differentiable vector A is: 


V- A = 


r 2 sin 9 


d(r 2 A r ) <9(sindA 0 ) dAg 

sm 9 y „ ’ + r — - — — — — + r— f 


dr 


d9 


(163) 


The curl of a differentiable vector A is: 


V x A = 


r 2 sin 9 


e r reg r sin 


_0 

dr 


d_ 

ae 


JL 

d(j> 


A r rAg r sin 9A ( j ) 


(164) 


For a general orthogonal system in a 3D space identified by the coordinates (ui,U 2 ,u 3 ) 


with unit basis vectors ui, 112 and U 3 and scale factors hi, hi and /13 where hi = 
r is the position vector: 


dr 


dui 


and 


1 ' Again, the summation convention does not apply to r, 0 and (j). 


4.4 Common Identities in Vector and Tensor Notation 


56 


The V operator is: 


ui d u-2 d u 3 d 

— 1 — ± 1 — ± 

h\ du\ h 2 du-2 h 3 du 3 


( 165 ) 


The Laplacian operator is: 


V 2 


1 

h\h 2 h 3 


d / h 2 h 3 d \ d_ ( hih 3 d \ _d_ ( h x h 2 d \ 

dui \ hi du\ ) du 2 V h 2 d u 2 ) du 3 \ h 3 du 3 ) 


( 166 ) 


The gradient of a differentiable scalar / is: 


V/ 


h i du i /i 2 dui h 3 du 3 


( 167 ) 


The divergence of a differentiable vector A is: 


V- A = 


hih2h 3 


r\ r\ r\ 

- — (h 2 h 3 Ai) + — — ( hih 3 A 2 ) + 7^ — (hih 2 A 3 ) 

OVL\ OXl^ ^^3 


( 168 ) 


The curl of a differentiable vector A is: 


V x A 


1 

hih 2 h 3 



h 2 u 2 

h 3 u 3 

d 

a 

a 

du\ 

du 2 

du 3 

hiAi 

/I2A2 

h 3 A 3 


( 169 ) 


4.4 Common Identities in Vector and Tensor Notation 


• Here we present some of the widely used identities of vector calculus in the traditional 
vector notation and in its equivalent tensor notation. In the following bullet points, / 
and h are differentiable scalar fields; A, B, C and D are differentiable vector fields; and 
r = Xie.i is the position vector. 
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V • r = n 

$ ( 170 ) 

d,x t = n 

where n is the space dimension. 

V x r = 0 

t ( 171 ) 

^■ijkdjXk 0 

V (a • r) = a 

$ ( 172 ) 

d i ( djXj ) = a* 

where a is a constant vector. 

V-(V/) = V 2 / 

$ ( 173 ) 

di ( dj ) = duf 
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• 

V ■ (V x A) = 0 

t 

^ijkdidj-A-k 0 

(174) 

• 

V x (V/) = 0 

t 

^ijkdjOkf 0 

(175) 

• 

VC/A) = fVh + hN7f 

t 

dt ( fh ) = fdih + hdj 

(176) 

• 

V • (/A) = /V-A + A-V/ 

$ 

A (/A) = ,/AA + AA/ 

(177) 
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V x (/A) = /VxA + V/xA 


$ 


( 178 ) 


^ ijkdj f ^ijkdj A-k “h €ijk -^-k 


A • (B x C) = C • (A x B) = B • (C x A) 

$ $ 


jC k ^kijC k-^-iB j €j ki B j C ^ 


( 179 ) 


Ax(BxC) = B (A ■ C) — C (A • B) 

t ( 180 ) 

tijkAjCklmBlCm B i (A 

m Cm) - Ci (ABt) 


Ax(VxB) = (VB) • A — A • VB 

t ( 181 ) 

£ijk£kim-AjdiB m = (9jB m ) A m Ai (diBi) 
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V x (V x A) = V (V • A) - V 2 A 

$ (182) 
( tj k ( k iw d)j ()/ .4 / , ) Dj (d m A m ) duAi 

V (A • B) = A x (V x B) + B x (V x A) + (A ■ V) B + (B ■ V) A 

t ( 183 ) 

(A m B m ) CjkAj (e k i m diB m ) + eij k Bj [ekimdiAm) H~ (A/c^) B j + (B/<9/) A j 

V • (A x B) = B • (V x A) - A • (V x B) 

$ (184) 

(cijkAjB k ) B k (e ki jdiAj) Aj (rji^OjBi.j 

Vx(AxB) = (B • V) A + (V • B) A — (V • A) B — (A • V) B 

t ( 185 ) 


Cijktkimdj ( AiB m ) — (B m d m ) Ai + (d m B m ) Ai — ( djAj ) Bi — ( Ajdj ) Bi 
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(A x B) • (C x D) 


AC AD 

B C B D 


t 


eijkAjBkeumCiDr 


( AiCi ) ( B m D m ) — ( A m D m ) ( BiCi ) 


(186) 


(AxB)x(CxD) = [D • (A x B)]C - [C • (A x B)]D 


t 


(187) 


tijk^jmnAmBnCkpqCpDq ( C qt nn B ( j /I rr I3 n ) ( ’ , y€pm n CpA m B n J Di 


• In vector and tensor notations, the condition for a vector field A to be solenoidal is: 


V-A = 0 


t 

diAi = 0 


(188) 


• In vector and tensor notations, the condition for a vector held A to be irrotational is: 


V x A = 0 




( 189 ) 


-Ak 


0 
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4.5 Integral Theorems in Tensor Notation 

• The divergence theorem for a differentiable vector field A in vector and tensor notation 
is: 



V • A dr = A n da 


v 


$ 


diAidr = / AiUidcr 


<v 


(190) 


where V is a bounded region in an nD space enclosed by a generalized surface S, dr and 
da are generalized volume and surface elements respectively, n and n l are unit normal to 
the surface and its i th component respectively, and the index % ranges over 1 , ,n. 

• The divergence theorem for a differentiable rank-2 tensor field A in tensor notation for 
the first index is given by: 

I diAudr = I Auriida (191) 

Jv J s 

• The divergence theorem for differentiable tensor fields of higher ranks A in tensor nota- 
tion for the index k is: 


9/x :Aij k m d>T / ...k...rnX^kdo' 


lv 


V- 


(192) 


Stokes theorem for a differentiable vector field A in vector and tensor notation is: 


(V x A) • nda = A • dr 

Jc 


$ 


(193) 


j eijkdjAkUida = j A t dxi 


where C stands for the perimeter of the surface S and dr is the vector element tangent to 
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the perimeter. 

• Stokes theorem for a differentiable rank-2 tensor held A in tensor notation for the first 
index is: 


^ '!] k ilj ' i kl ( ( 1 drT 


Audxi 


(194) 


Js Jc 

• Stokes theorem for differentiable tensor fields of higher ranks A in tensor notation for 
the index k is: 






n,do = A 




dxi 


(195) 


'C 


4.6 Examples of Using Tensor Techniques to Prove Identities 


V ■ r = diXi 

(Eq. 139) 

= Su 

(Eq. 80) 

= n 

(Eq. 80) 

[V x r] . = ei jk djX k 

(Eq. 146) 

^ijk^kj 

(Eq. 79) 

= e ijj 

(Eq. 76) 

= 0 

(Eq. 71) 


(196) 


( 197 ) 


Since i is a free index the identity is proved for all components. 
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• V (a • r) = a: 


[V (a ■ r)] . = di ( djXj ) 

(Eqs. 134 & 120) 

= djdiXj + XjdjUj 

(product rule) 

= CLjdiXj 

(a j is constant) 

= a,j5ji 

(Eq. 79) 

CL i 

(Eq. 76) 

= [a]i 

(definition of index) 

Since % is a free index the identity is proved for all components. 
. V • (V/) = V 2 /: 

V ■ (V/) = d , f [V/] 4 

(Eq. 139) 

= di ( dif ) 

(Eq. 134) 

= did if 

(rules of differentiation) 

= duf 

(definition of 2nd derivative) 

= v 2 / 

(Eq. 151) 


(198) 


( 199 ) 
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• V • (V x A) = 0: 

V • (V x A) = di [V x A] . 

di (gi.jkdj yfk) 
^ijkdidj A k 

^ijkdjdiA k 

EjikdjdiA k 

€ijkdidjA k 

= 0 


(Eq. 139) 

(Eq. 146) 

(i d not acting on e) 

(continuity condition) (200) 

(Eq. 84) 

(relabeling dummy indices i and j) 

(since €ij k d{dj A k e %j k di dj A k ) 


This can also be concluded from line three by arguing that: since by the continuity con- 
dition di and dj can change their order with no change in the value of the term while a 
corresponding change of the order of i and j in e l]k results in a sign change, we see that 
each term in the sum has its own negative and hence the terms add up to zero (see Eq. 
94). 

. V x (V/) = 0: 

[V X (V/)], = e ijk d 3 [V/] fe (Eq. 146) 

= tijkdj ( d k f ) (Eq. 134) 


= €ij k djdkf (rules of differentiation) 

= €ij k d k djf (continuity condition) 


( 201 ) 


= -Ukjdkdjf (Eq. 84) 


^ ijkdjd k f 


(relabeling dummy indices j and k ) 


= 0 


(since e ijk djd k j = - e ijk djd k f ) 


This can also be concluded from line three by a similar argument to the one given in the 
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previous point. Because [V x (V/)]j is an arbitrary component, then each component is 
zero. 

•V(fh) = fVh + hVf: 


[V (/A)], = a, ( fh ) 

(Eq. 134) 

= fdji + hdif 

(product rule) 

(202) 

= [/VfcL + [hVf], 

(Eq. 134) 

= [/Vfe + ftV/] ( 

(Eq. 16) 


Because i is a free index the identity is proved for all components 

. V-(/A) = /V-A + A- V/: 


V ■ (/A) = ^ [/A] i 

(Eq. 139) 

= di ( fAi ) 

(definition of index) 

(203) 

= fdiAi + Aid if 

(product rule) 

= /V ■ A + A ■ V/ 

(Eqs. 139 & 153) 


• V x (/A) = /V x A + V/ x A 


[v x (/A)], = [/A], 

(Eq. 146) 

^ijk^j (f Ah') 

(definition of index) 

f ^ijk^j Ak + €-ijk ifijf') A-k 

(product rule & commutativity) 

(204) 

f ^ijk^j A-k 4“ ^-ijk [ ^ /] j -A-h 

(Eq. 134) 

= [/V x A], + [V/ x A], 

(Eqs. 146 & 121) 

= [/V x A + V/ x A], 

(Eq. 16) 


Because i is a free index the identity is proved for all components. 
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• A • (B x C) = C • (A x B) = B • (C x A): 


Bj C k 

(Eq. 122) 

£ kij A{ Bj C k 

(Eq. 84) 

tkijCfcAi Bj 

(commutativity) 

C • (A X B) 

(Eq. 122) 

€j ki-Ai Bj C k 

(Eq. 84) 

ki Bj CkAi 

(commutativity) 

B • (C x A) 

(Eq. 122) 


(205) 


The negative permutations of these identities can be similarly obtained and proved by 
changing the order of the vectors in the cross products which results in a sign change. 
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• A x (B x C) = B (A ■ C) - C (A • B): 


[A x (B x C)] 4 = e ijk Aj [B x C] fc 

(Eq. 121) 

^ ijk Ay ( k/ jn P ! ( 'm 

(Eq. 121) 

^ijk^klmAj BiC m 

(commutativity) 

( ijk ( l tn k - P ! ( 'in 

(Eq. 84) 

( jl ) AjBiC m 

(Eq. 102) 

fiiifijmAjBiC m 8 im 8jiAjBiC m 

(distributivity) 

(SilBi) (8j m AjC m ) (8 irr/C rn . ) ( 8jl Aj Bj ) 

(commutativity and grouping) 

= Bi (. A m C m ) — Ci ( AiBi ) 

(Eq. 76) 

= Bi (A • C) - Q (A • B) 

(Eq. 120) 

= [B (A • C)] ■ — [C (A • B)] . 

(definition of index) 

— [B (A • C) — C (A • B)] i 

(Eq. 16) 

(206) 


Because i is a free index the identity is proved for all components. Other variants of this 
identity [e.g. (A x B) x C] can be obtained and proved similarly by changing the order 
of the factors in the external cross product with adding a minus sign. 
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A x (V x B) = (VB) A - A VB 


[A x (V x B)]. = e ijk Aj [V x B] k 

(Eq. 121) 

£ijkAj6klmdlBm 

(Eq. 146) 

^ ij k t kirn Aj Oj B rn 

(commutativity) 

^ijk^lmkAjOlB m 

(Eq. 84) 

AjdiB m 

(Eq. 102) 

SiiSjmAjdiB m 5im5jiAjdiB m 

(distributivity) 

= A m diB m - AidiBi 

(Eq. 76) 

= ( diB m ) A m - Ai ( d t Bi ) 

(commutativity & grouping) 

= [(VB) • A] . - [A • VB] i 

(Eq. 135 & § 2.6.5) 

= [(VB) • A — A • VB] . 

(Eq. 16) 

(207) 


Because i is a free index the identity is proved for all components. 
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• V x (V x A) = V (V • A) - V 2 A: 


[V x (V x A)] ( = [V x A]„ 

(Eq. 146) 

£ ijkdj \^klm9lA m j 

(Eq. 146) 


(' d not acting on e) 

^ijk^lmkdjdiA m 

(Eq. 84 & definition of derivative) 

9jdiA m 

(Eq. 102) 

jmdjdiA rn 5i m SjidjdiA rn 

(distributivity) 

= dmdiAm - didlAi 

(Eq. 76) 

= di ( dmAm ) - dllAi 

(' d shift, grouping & Eq. 4) 

— [V (V • A)]^ — [V 2 A] . 

(Eqs. 139, 134 & 152) 

= [V (V • A) — V"A] i 

(Eqs. 16) 


(208) 


Because i is a free index the identity is proved for all components. This identity can also 
be considered as an instance of the identity before the last one, observing that in the 
second term on the right hand side the Laplacian should precede the vector, and hence no 
independent proof is required. 

• V (A • B) = A x (V x B) + B x (V x A) + (A • V) B + (B ■ V) A: 
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We start from the right hand side and end with the left hand side 

[A x (V x B) + B x (V x A) + (A ■ V) B + (B • V) A] i 
[A x (V x B)] i + [B x (V x A)], + [(A • V) B] i + [(B • V) A] ■ 
e ijk-Aj [V x B]j. + f-ijkBj [V x A] fe + (A;d;) Bi A ( Bidi ) A, 
Uhl m diB m ) + e ijkBj ( €kl m 8 l A m ) + (A l d l )B i + (B l d l )A l 
€ ijk^klmAjdiB rn A 6 ijk^-kl r nBjdiA rn A ( A i ()j ) Bi A ( Bj 0 [ j A 7 
tijk^lmkAjOlBm A ^ijkV-lmkB jd[A rn A ( A / Oj j Bj A ( B/ Dj ) A t 
(fiilfijm firrnAjl ) Aj ()j B rn A {SilSjm BjdlA m A ( A / (); ) Bj A ( Bj D ; j Aj 

{,^ii^jvaAjdiB m t) vm Sji Aj ()i B rn ) A (SuSjmBjdiAn A l/rn Sji Bj [)/ A rn j A ( A jDj) Bj A ( Bj Oi ) Aj 
SudjmAjdiB m — AjdjBj A SnSjmBjdiAm — BidiAj A ( Aidi ) Bj A ( Bidi ) Aj 
SuSjmAjdiBm — ( Aidi ) Bj A SjiSj m BjdiA m — ( Bjdi ) Aj A ( Aidi ) Bj A ( Bjdi ) Aj 

SuSjmAjdiBn A 6 u 6 jm BjdiA m 
A m d, B rn A B m djA m 
dj ( A m B m ) 


(Eq. 16) 

(Eqs. 121, 139 & indexing) 
(Eq. 146) 

( commut at i vity ) 

(Eq. 84) 

(Eq. 102) (209) 

(distributivity) 

(Eq. 76) 

(grouping) 

(cancellation) 

(Eq. 76) 

(product rule) 

[V (A • B)] j (Eqs. 134 & 139) 


Because i is a free index the identity is proved for all components. 
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• V • (A x B) = B • (V x A) - A • (V x B): 


V • (A x B) = di [A x B] . 

9i ( UjkAj 13 k ) 

C ijkdi (^AjBk) 

£ijk ( B kd/ Aj -f- AjdiBk) 

^ ij h B} { ()j Aj -\~ k. jj k Aj Oj B [. 

^kij B k d,Aj e Tl kAjdjB k 

Bk ( tkijdiyf j ) Aj (ejikdiBk) 

= B k [V x A] k - Aj [V x B]^ 

= B ■ (V x A) — A ■ (V x B) 


(Eq. 139) 

(Eq. 121) 

(d not acting on e) 

(product rule) 

(distributivity) (210) 

(Eq. 84) 

(commutativity & grouping) 


(Eq. 146) 
(Eq. 120) 
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• V x (A x B) = (B ■ V) A + (V • B) A - (V • A) B - (A ■ V) B: 

[V x (A x B)] i = e ijk dj [A x B} k 

C ijkdj 

^ijk^klmdj 

tijk^klm (B m dj A/ “I - AidjB m ) 
tijk^lmk {B m djAi A AidjB m ) 

i&il&jm ( B m Oj A / A Aidj -B m ) 

jm.Bmff j Ai A hjj d j fn Aj Oj B rn Si m SjiB m djAi 5 im 5jiAidjB m 

— B m d m Ai A Aid m B m — BidjAj — AjdjBi 


(Eq. 146) 

(Eq. 121) 

(' d not acting on e 
(product rule) 

(Eq. 84) 

(Eq. 102) 
(distributivity) 
(Eq. 76) 


— (. B m d m ) Ai A (d m B m ) Ai — ( djAj ) Bi — ( Ajdj ) Bi 
= [(B • V) A], A [(V • B) A] i - [(V ■ A) B] . - [(A • V) B] . 

= [(B • V) A A (V ■ B) A - (V • A) B - (A • V) B] . 


(grouping) 

(Eqs. 153 & 139) 


Because i is a free index the identity is proved for all components. 


(Eq. 16) 
( 211 ) 
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• (A x B) ■ (C x D) 


AC AD 
B C B D 


(AxB).(CxD) = [AxB] j [CxD] j 

tijkAjBkEiimCiD 

m 

^■ijk^-ilmAjB} t ClE rn 

if j A km fjrn.fkl ) r\ j Up. Cl] E m 

fjlfikm Aj Pk f f E rn fjm Sj : } Aj 13 j. Cj E rn 

if l ) i/i kn/P kErn) if jm A j E rn j i fkl E k.Ul ) 

= (A t Ci) ( B m E m ) - (A m D m ) (BtQ) 

= (A ■ C) (B • D) — (A • D) (B-C) 

AC AD 
B C B D 


(Eq. 120) 

(Eq. 121) 

(commutativity) 

(Eqs. 84 & 102) 
(distributivity) 
(commutativity & grouping) 
(Eq. 76) 

(Eq. 120) 

(definition of determinant) 
( 212 ) 
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. (A x B) 

x(CxD) 

= [D • (A x B)] C - [C • (A x B)] D: 


[(A x B) 

x (C x D)] ( 

= e ijk [A x B] j [C x D] fc 

(Eq. 121) 



£ijk€jmn-ArnB n €kpqCpDq 

(Eq. 121) 



£ijk£kpq£jmn-A-mB n CpDq 

(commutativity) 



tijkCpqktjmn-AmBjiCpDq 

(Eq. 84) 



ifiip^jq $iq$jp) Cjmn-A-mBnCpDq 

(Eq. 102) 



(fiipfijqtjmn ^iq^jp^jmn) A mB n CpDq 

(distributivity) 



y^ip^qmn ^iq^-pmn) A mB n CpD q 

(Eq. 76) 



$ip€qmnA m B n CpD q 5iq€p Tnn A rn B n CpD q 

(distributivity) 



£qmn-A-rnB n C iD q 6p Tnn A rn B n CpD-i 

(Eq. 76) 



^qrnnB qA^B n C i €p mn C pA m B n D ^ 

(commutativity) 



\€qmnB qA m B n j C{ D} 

(grouping) 



= [D • (A x B)] Ci - [C • (A x B)] A 

(Eq. 122) 



= [[D • (A x B)] C]. - [[C • (A x B)] D]. 

(definition of index) 



— [[D • (A x B)] C — [C • (A x B)] D] . 

(Eq. 16) 


( 213 ) 


Because i is a free index the identity is proved for all components 
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5 Metric Tensor 

• This is a rank-2 tensor which may also be called the fundamental tensor. 

• The main purpose of the metric tensor is to generalize the concept of distance to gen- 
eral curvilinear coordinate frames and maintain the invariance of distance in different 
coordinate systems. 

• In orthonormal Cartesian coordinate systems the distance element squared, ( ds ) 2 , be- 
tween two infinitesimally neighboring points in space, one with coordinates x l and the 
other with coordinates x l + dx l , is given by 

(ds) 2 = dx l dx l = 8 i jdx l dx :i (214) 

This definition of distance is the key to introducing a rank-2 tensor, g t j , called the metric 
tensor which, for a general coordinate system, is defined by 

(ds) 2 = gijdx l dx J (215) 

The metric tensor has also a contravariant form, i.e. g l T 

• The components of the metric tensor are given by: 

gij = E i • E j & g ij = E i ■ E j (216) 

where the indexed E are the covariant and contravariant basis vectors as defined in § 2.5.1. 

• The mixed type metric tensor is given by: 

g) = E ■ Ej = S) & g/ = E • E = 5/ (217) 


and hence it is the same as the unity tensor. 


5 METRIC TENSOR 


77 


• For a coordinate system in which the metric tensor can be cast in a diagonal form where 
the diagonal elements are ±1 the metric is called flat. 

• For Cartesian coordinate systems, which are orthonormal flat-space systems, we have 

9 ij = S ij = 9ij = 5ij (218) 

• The metric tensor is symmetric, that is 

9 ij = 9 ji & g' 3 = 9 " (219) 

• The contravariant metric tensor is used for raising indices of covariant tensors and the 

covariant metric tensor is used for lowering indices of contravariant tensors, e.g. 

A { = g ij Aj At = gijA j (220) 

where the metric tensor acts, like a Kronecker delta, as an index replacement operator. 
Hence, any tensor can be cast into a covariant or a contravariant form, as well as a mixed 
form. However, the order of the indices should be respected in this process, e.g. 

AT = g jk A ik ± A? = g jk A ki (221) 

Some authors insert dots (e.g. A’?) to remove any ambiguity about the order of the indices. 

• The covariant and contravariant metric tensors are inverses of each other, that is 

[9«]=[s y ] _1 & [S ii ]=\gn r‘ ( 222 ) 

Hence 

9 ik 9 kj = S ) & g ik g ki S, 1 (223) 
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• It is common to reserve the “metric tensor” to the covariant form and call the con- 
travariant form, which is its inverse, the “associate” or “conjugate” or “reciprocal” metric 
tensor. 

• As a tensor, the metric has a significance regardless of any coordinate system although 
it requires a coordinate system to be represented in a specific form. 

• For orthogonal coordinate systems the metric tensor is diagonal, i.e. g tJ = g = 0 for 

• ^ j- 

• For flat-space orthonormal Cartesian coordinate systems in a 3D space, the metric tensor 
is given by: 


[dij] [$ij] 


1 0 0 
0 1 0 
0 0 1 


~5 ij ] = [g ij ' 


(224) 


• For cylindrical coordinate systems with coordinates (p,(f>,z), the metric tensor is given 
by: 


1 

0 

0 



1 

0 

0 

0 

P 2 

0 

& 

W = 

0 

J_ 

p 2 

0 

0 

0 

1 



0 

0 

1 


(225) 


• For spherical coordinate systems with coordinates (r, 6 , (f ) ), the metric tensor is given by: 


1 

0 

0 



l 

0 

0 

0 

r 2 

0 

& 

m = 

0 

1 

r 2 

0 

0 

0 

r 2 sin 2 6 



0 

0 

1 

r 2 sin 2 6 


( 226 ) 
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6 Covariant Differentiation 


• The ordinary derivative of a tensor is not a tensor in general. The objective of covariant 
differentiation is to ensure the invariance of derivative (i.e. being a tensor) in general 
coordinate systems, and this results in applying more sophisticated rules using Christof- 
fel symbols where different differentiation rules for covariant and contravariant indices 
apply. The resulting covariant derivative is a tensor which is one rank higher than the 
differentiated tensor. 

• Christoffel symbol of the second kind is defined by: 


(k \ _ g kl f dgu d 9jl dg,A 

2 y dx ] dx l dx l ) 


(227) 


where the indexed g is the metric tensor in its contravariant and covariant forms with 
implied summation over /. It is noteworthy that Christoffel symbols are not tensors. 

• The Christoffel symbols of the second kind are symmetric in their two lower indices: 


ttl = {%) 


(228) 


• For Cartesian coordinate systems, the Christoffel symbols are zero for all the values of 
indices. 

• For cylindrical coordinate systems ( p,4>,z ), the Christoffel symbols are zero for all the 
values of indices except: 


{y = -P (229) 

{£.} = {!.} = l 


where (1,2,3) stand for (p,(f>,z). 
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• For spherical coordinate systems (r,0,(f>), the Christoffel symbols are zero for all the 
values of indices except: 


ft} = 

— r 


{y = 

—r sin 2 9 


(y = 

ft} = 

1 

r 

&} = 

— sin 9 cos 9 

(u = 

ft} = 

1 

r 

<y = 

ft} = 

cot 


(230) 


where (1,2,3) stand for ( r,9,(j ) ). 

• For a differentiable scalar / the covariant derivative is the same as the normal partial 
derivative, that is: 

U = f, = dj (231) 


This is justified by the fact that the covariant derivative is different from the normal partial 
derivative because the basis vectors in general coordinate systems are dependent on their 
spatial position, and since a scalar is independent of the basis vectors the covariant and 
partial derivatives are identical. 

• For a differentiable vector A the covariant derivative is: 



SiA, ~ {‘ } A 

3i V + {!.} A 


k 


(covariant) 

(contravariant) 


( 232 ) 
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• For a differentiable rank- 2 tensor A the covariant derivative is: 


Ajk-,i { j i} Aik {fci} Aji 

(covariant) 


A% = SiA* + {'} A lk + {‘} A* 

(contravariant) 

(233) 

A% = d,A] + {*} A\ - {',} A* 

(mixed) 



• For a differentiable rank-n tensor A the covariant derivative is: 

Aij.-.k o A ij-k I (i \ Adj...k , (j \ Aia...k _i_ . . . i fk \ Aij-a tOQ/L'l 

■ tS ~lm...p,q u q- ti 'lm...p ' \ aq I ' \ aq I ^ ' \ aql ■ n -lm...p 

_/«l Aij—k _ fa \ Aij-k _ . . . _ fa \ Aij-k 
\lqf' rl am...p \mq / ' rl la...p l pq) 

• From the last three points a pattern for covariant differentiation emerges, that is it starts 
with a partial derivative term then for each tensor index an extra Christoffel symbol term 
is added, positive for superscripts and negative for subscripts, where the differentiation 
index is the second of the lower indices in the Christoffel symbol. 

• Since the Christoffel symbols are identically zero in Cartesian coordinate systems, the 
covariant derivative is the same as the normal partial derivative for all tensor ranks. 

• The covariant derivative of the metric tensor is zero in all coordinate systems. 

• Several rules of normal differentiation similarly apply to covariant differentiation. For 
example, covariant differentiation is a linear operation with respect to algebraic sums of 
tensor terms: 

d-i (aA ± 6B) = ad-iA ± bdjD (235) 

where a and b are scalar constants and A and B are differentiable tensor fields. The 
product rule of normal differentiation also applies to covariant differentiation of tensor 
multiplication: 


d,i (AB) = (d ;i A) B + A<9 ;i B 


(236) 
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This rule is also valid for the inner product of tensors because the inner product is an 
outer product operation followed by a contraction of indices, and covariant differentiation 
and contraction of indices commute. 

• The covariant derivative operator can bypass the raising/lowering index operator: 

A = fhjA 3 => d- m Ai = (i l] d ]m A 3 (237) 

and hence the metric behaves like a constant with respect to the covariant operator. 

• A principal difference between normal partial differentiation and covariant differentiation 
is that for successive differential operations the partial derivative operators do commute 
with each other (assuming certain continuity conditions) but the covariant operators do 
not commute, that is 


d idj = djdi but d-id.j ^ d :J d :l (238) 

• Higher order covariant derivatives are similarly defined as derivatives of derivatives; 
however the order of differentiation should be respected (refer to the previous point). 



REFERENCES 


83 


References 

[1] G.B. Arfken; H. J. Weber; F.E. Harris. Mathematical Methods for Physicists A Com- 
prehensive Guide. Elsevier Academic Press, seventh edition, 2013. 

[2] R.B. Bird; R.C. Armstrong; O. Hassager. Dynamics of Polymeric Liquids , volume 1. 
John Wiley & Sons, second edition, 1987. 

[3] R.B. Bird; W.E. Stewart; E.N. Lightfoot. Transport Phenomena. John Wiley & Sons, 
second edition, 2002. 

[4] M.L. Boas. Mathematical Methods in the Physical Sciences. John Wiley & Sons Inc., 
third edition, 2006. 

[5] C.F. Chan Man Fong; D. De Kee; P.N. Kaloni. Advanced Mathematics for Engineering 
and Science. World Scientific Publishing Co. Pte. Ltd., first edition, 2003. 

[6] T.L. Chow. Mathematical Methods for Physicists: A concise introduction. Cambridge 
LIniversity Press, first edition, 2003. 

[7] J.H. Heinbockel. Introduction to Tensor Calculus and Continuum Mechanics. 1996. 

[8] D.C. Kay. Schaum’s Outline of Theory and Problems of Tensor Calculus. McGraw- 
Hill, first edition, 1988. 

[9] K.F. Riley; M.P. Hobson; S.J. Bence. Mathematical Methods for Physics and Engi- 
neering. Cambridge LIniversity Press, third edition, 2006. 

[10] D. Zwillinger, editor. CRC Standard Mathematical Tables and Formulae. CRC Press, 
32nd edition, 2012. 

Note: As well as the references cited above, I benefited during the writing of these notes 
from many sources such as tutorials, presentations, and articles which I found on the 



REFERENCES 


84 


Internet authored or composed by other people. As it is difficult or impossible to retrace 
and state all these sources, I make a general acknowledgment to all those who made their 
documents available to the public. 



